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We investigate pairing symmetry and transition temperature in a quasi-one-dimensional re- 
pulsive Hubbard model. We solve the Eliashberg equation using the third-order perturbation 
expansion with respect to the on-site repulsion U . We find that when the electron number den- 
sity is shifted from the half-filled state, a transition into unconventional superconductivity is 
expected. When one-dimensionality is weak, a spin-singlet state is favorable. In contrast, when 
one- dimensionality is strong and electron number density is far from the half-filled state, a spin- 
triplet state is stabilized. Finally, we discuss the possibility of unconventional superconductivity 
caused by the on-site Coulomb repulsion in /3-Nao.33V205. 

KEYWORDS: quasi-one-dimensional superconductors, pairing symmetry, transition temperature, third- 
order perturbation theory 



Superconductivity in quasi-one-dimensional conduc- 
tors has been studied as an important phenomenon. To- 
day, some quasi-one-dimensional superconductors, such 
as (TMTSF)2Xi^2 ^^^^ Sri4_:ECa:,Cu2404i,^ have been 
discovered, and their superconductivity has been inves- 
tigated. Recently, a superconducting transition in /3- 
Nao.33V205, which has a quasi-one-dimensional lattice 
structure similar to Sri4_2:Caa;Cu2404i, has been dis- 
covered,^ and this phenomenon has attracted our atten- 
tion. The transition temperature is Tc ~ 8 K under high 
pressures of approximately 8 GPa. At ambient pressure, 
this material shows quasi-one-dimensional metallic be- 
havior in an electric resistivity experiment at high tem- 
perature,^ and encounters a charge-ordered transition at 
Tco — 135 K.^ Furthermore, below Tn ~ 25K, a antifer- 
romagnetic ordered phase appears in the charge-ordered 
phase. ^ Under high pressures of approximately 8 GPa, 
the charge-ordered phase abruptly vanishes, accompa- 
nied by the superconducting transition.'* It is not clear 
under high pressure whether the antiferromagnetic phase 
in the charge-ordered phase survives or not. However, the 
existence of the antiferromagnetic phase at ambient pres- 
sure suggests that the electron correlation is important. 

Such an electron correlation effect leads to uncon- 
ventional superconductivity rather than conventional s- 
wave superconductivity induced by the electron-phonon 
coupling.^ Such investigations have already been re- 
ported in the quasi-one-dimensional superconductors, 
(TMTSF)2X and Sri4_j;Caa:Cu2404i. The superconduc- 
tivity in (TMTSF)2X has been investigated using the 
fluctuation-exchange approximation (FLEX)^ and the 
third-order perturbation theory (TOPT).^ Both theoret- 
ical calculations suggest that a d-wave like spin-singlet 
state is the most stable. Also, in Sri4_a;Caj;Cu2404i, 
the FLEX calculation for the trellis lattice indicated 
a d-wave like spin-singlet state.*" However, experimen- 
tally, in both materials, the Knight shift does not change 
above and below Tc, and the spin-triplet state is indi- 
cated,**'*^ although it is still confusing. Thus, intensive 
investigations on quasi-one-dimensional superconductors 
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Fig. 1. (a) Schematic figure of tlie lattice used in this calcula- 
tion. < « < 7) is the hopping integral. The region enclosed 
by rectangle is primitive cell. The primitive cell topologycally 
composes triangular lattice, (b) The Fermi surface for n = 0.90, 
to = 1.0. Since the lattice is topologycally triangular lattice, the 
Brillouin zone is hexagonal. 



should be carried out. In this letter, we investigate in 
detail unconventional superconductivity in a quasi-one- 
dimensional Hubbard model, taking up the superconduc- 
tivity in /3-Nao.33V205. 

Now, let us consider the lattice structure and the band 
structure in /3-Nao.33V205. This material has three types 
of vanadium site, VI, V2 and V3. VI is on the VOe 
zigzag chain, V2 is on the VOe ladder chain and V3 is 
on the VO5 zigzag chain. At ambient pressure, NMR 
experiments indicated that V3 does not have any con- 
duction electrons at low temperatures.*'^ Assuming that 
conduction electrons at V3 are also empty under high 
pressures, we can consider the lattice structure with the 
V network shown in Fig. 1(a), which is important for 
electric conductivity. This structure is different from the 
trellis lattice only in the VI zigzag chain in the middle 
of Fig. 1(a). The unit cell is a thick- line rectangle, and 
contains four V sites. Since there is no information on 
the band structure of /3-Nao.33V205, we discuss a sim- 
ple tight-binding model with an s orbital on each V site. 
In this case, we consider 7 types of hopping integrals 
(^1 ^ tr) displayed in Fig. 1(a). By numerically diago- 
nalizing the energy matrix, we obtain four bands, since 
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there arc four orbitals in the unit cell. Among them, wc 
only use the lowest energy band. In /3-Nao.33V205, there 
is one electron per unit cell, if we simply count the va- 
lence electrons. Under the ideal condition that all elec- 
trons occupy only the lowest energy band, it becomes 
half-filled. Therefore, we deal with the electron number 
density n as a parameter less than the half-filled state. If 
the superconductivity of this material is caused by elec- 
tron correlation, electron correlation must be strong. For 
electron correlation to be strong, a high-density band is 
required. Therefore, the electrons should mainly occupy 
the lowest energy band. Prom these viewpoints, we use 
the single-band model. If there is no high-density band, 
we have to consider a mechanism other than electron cor- 
relation. As a typical set of parameters, we use ti = 1.0, 
t2 = to, ts = t4 = t5 = Q.Sto and = = 0.2to- Here, to 
is a measure of one-dimensionality, and with decreasing 
to, the Fermi surface becomes more one-dimensional. We 
assume that the results of calculations mainly depend 
on the one-dimensionality of the Fermi surface. Actu- 
ally, when we change the ratio of the transfer integrals 
ti maintaing the one-dimensionality of the Fermi sur- 
face, the results of the calculations are almost unchanged. 
Therefore, we can study the dependence of the form of 
the Fermi surface using the parameter to. In Fig. 1(b), 
we show a typical quasi-one-dimensional Fermi surface 
for to = 1.0 and the electron number density n = 0.90. 
This Fermi surface possesses a less nesting property, and 
is different from the band structures with an almost per- 
fect nesting property discussed so far in the quasi-one- 
dimensional model calculations. Here, we investigate in 
detail superconductivity in such a situation. We consider 
the quasi-one-dimensional single-band Hubbard model 
with the lowest band e{k) discussed above, 
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We treat this model using the third-order perturbation 
expansion. Hereafter, in order to obtain a moderate tran- 
sition temperature Tc, wc set IJ = 5.0, which is almost 
equal to the bandwidth. The third-order perturbation 
theory in the strongly correlated region has been justi- 
fied by higher-order calculations of pairing interactions."'^^ 
We can apply the perturbation theory for the appropri- 
ate values of U to obtain the reliable value of T^. 

We apply the third-order perturbation theory with re- 
spect to U to our model. Diagrams of the normal self- 
energy are shown in Fig. 2. The normal self-energy is 
given by 



SN(fc) = - Yy\o{h - k')Go{k') 



(2) 
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Fig. 2. The diagrams of Normal self-energy. The solid lines rep- 
resent Green's function Go(fe) and the broken lines represent the 
Coulomb repulsion U. 



where 



Go(fc) 



iuJn - s{k) + IJ,' 



(«) = -iv5I^o(^)^o(g + fc), 



(3) 



'o{q) = -j;^Y.(^o{k)Goiq-k). 



Here, Go{k) with the short notation k = {k,uin) rep- 
resents the bare Green's function. Since the first-order 

normal self-energy is constant, it can be included by the 
chemical potential jj.. The dressed Green's function G{k) 
is given by 



G{k) 



(4) 



iu>n - e{k) - I]N(fc) + /'' + (^/i ' 

Here, the chemical potential and the chemical potential 
shift (5/x are determined so as to fix the electron number 
density n, 

^ = ^^T.Go{k) = 2^J2Gik)- (5) 

k k 

We also expand the effective pairing interaction up to 
the third order with respect to U. For the spin-singlet 
state, the effective pairing interaction is given by 

^Singiet(^. fc,) ^ Vi^f'\k; k') + V^:Xt\k; k'), (6) 



where 

Vl^t\k; k') = U + U\o{k - k') + 2U\l{k - k'), 



(7) 



and 



Vv^'^ik;k') = 2{T/N)Re[Y,Go{ki) 



(8) 

X {xo{k + ki) - Mk + ki))Go{k + fci - k')U^' . 
For the spin-triplet state, 

y™piet(fc. fc,) = V^^t\k; k') + V^^^l^{k; k'), (9) 
where 

V^^t\k;k') = -U\o{k-k'), (10) 
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Fig. 3. Calculated maximum eigenvalues Amax for spin-singlet (or 
spin-triplet) state. The line with white squares (circles) is the 
result for the spin-singlet (spin-triplet) state obtained using the 
third-order perturbation theory. The line with the black squares 

(circles) is the result for the spin-singlct (spin-triplet) state with- 
out the pairing interaction due to the third-order terms. The 
parameters are n = 0.90 and to = 1.0. 



and 



(11) 

X (xo(fc + fci) + + ki))Go{k + fci - fc')f^^^ 



Here, K 
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(fc, k') is called the RPA terms and 

(fc, k') is called the vertex corrections. 
Near the transition point, the anomalous self-energy 
A(fc) satisfies the linearized Eliashberg equation, 

A^axA(fc) = E k')\G{k')\^^{k'), (12) 

k' 

where, V{k]k) is l/Singiet^-^. ^./^ yTripict(-^,. ^/^^ 
Amax is the largest positive eigenvalue. Then, the tem- 
perature at Amax = 1 corresponds to T^. By estimating 
Amax: wc Can determine which type of pairing symmetry 
is stable. For numerical calculations, we take 128 x 128 
fe-meshes for twice space of the first Brillouin zone and 
2048 Matsubara frequencies. 

In Fig. 3, we show the results for Amax in the case 
with n = 0.90 and to = 1.0. With decreasing tempera- 
ture, Amax increases. The spin-singlet state and the spin- 
triplet state possess almost the same transition tempera- 
ture, Tc ~ 0.004. If we assume that the bandwidth ~ 5 
corresponds to 1 eV, then ^ 8 K is obtained in accor- 
dance with the experimental value for /3-Nao.33V205. In 
Fig. 3, we also show the results for Amax obtained with- 
out the pairing interaction due to the third-order terms. 
For the spin-triplet state, we see that the vertex correc- 
tions are important for stabilizing the spin-triplet state 
from the comparison. 

In Fig. 4, wc show the momentum dependence of 
Xo(Q) 0) and </'o('7, 0) in the case with n = 0.90, to = 1.0 
and T = 0.004. Since the Fermi surface has a nest- 
ing property, Xo(q,0) has peaks near q = {tt,tt) and 
q = (0, 7r). Since the Brillouin zone is hexagonal, (7r,7r) 
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Fig. 4. (a) Contour plots of Xo(9iO). The peaks exist near q = 

(7r,7r) and q = (0, tt) (b) Contour plots of (f)o{q,0). The peaks 
exist at q = (0, 0). The parameters are n = 0.90, to = 1.0 and 
T = 0.004. 
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Fig. 5. (a) Contour plot of the anomalous self-energy. A(fc) = 
for the spin-singlet state has d-wave like momentum dependence, 
(b) Contour plot of the anomalous self-energy. A(fc) = for the 
spin-triplet state has p-wave like momentum dependence. The 
parameters are n = 0.90, to = 1.0 and T = 0.01. 



is equivalent to (0, tt). At the half-filled state, Xo{q,0) 
has peaks beside q = (tTjTt) and q = (0,7r). If the elec- 
tron number density n is shifted from the half-filled state, 
then the peaks are shifted from q= (tt, tt) and q = (0, tt) 
. On the other hand, 4'o{q, 0) has a peak a.t q = (0, 0) . 

In Fig. 5, we show the contour plots of the anomalous 
self-energy in the case of n = 0.90, to = 1.0 and T = 0.01. 
For the spin-singlet state, the momentum dependence of 
the anomalous self-energy on the Fermi surface is a d- 
wave like state with node. For the spin-triplet state, the 
momentum dependence of the anomalous self-energy is 
p-wave like, and is a fully gapped state. For the spin- 
singlet state, the RPA terms are dominant in the efii'ective 
interaction terms. In this case, we can easily understand 
the gap structure in Fig. 5(a) from the structure in the 
Eliashberg equation as follows. The peak structures of 
Xo(q, O) in Fig. 4(a) originate from the nesting property, 
between around point A and around point B, or around 
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Fig. 6. Amax as a function of n. The parameters are T = 0.01 and 
to = 1.0. 
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Fig. 7. Amax as a function of to. The parameters are n = 0.90 
and T = 0.01. 

point B' on the Fermi surfaces in Fig. 5(a). In order to 

obtain a positive value of Anjaxj it is favorable that signs 
of A (A;) around points B and B' are different from its 
sign around point A. The structure of A(fc) in Fig. 5(a) 
just becomes so. 

In Fig. 6, n dependence of Amax is shown. If n is 
near the half-filled state, the spin-singlet and the spin- 
triplet states have low values of Amax- With decreasing 
n, Amax increases, and the spin-singlet and spin-triplet 
states yield nearly the same Amax at n ~ 0.90. More- 
over, with decreasing n from n = 0.90, the values of 
Amax for the spin-triplet state become larger than those 
for the spin-singlet state. This indicates that the spin- 
triplet state may be realized far from the half-filled state. 
We can easily understand why the values of Amax are 
suppressed for the spin-singlet and spin-triplet states at 
around the half-filled state. For spin-singlet state, at the 
half-filled state, Xo(9jO) has a peak beside q — (0,7r). 
Considering the structure of the Eliashbcrg equation, 
in order to obtain a large positive value of Amax, it is 
not favorable that signs of A(fc) around point A are the 
same as its sign around point C on the Fermi surface 
in Fig. 5(a). Therefore, the values of Amax are strongly 
suppressed around the half-filled state by the conflicting 



peaks of Xo{q, 0). On the other hand, if the electron num- 
ber density is far from the half-filled state, the peak of 
Xo{q, 0) is far from q = (0, n). Therefore, the suppression 
of the values of Amax becomes weak. For the spin-triplet 
state, when the Fermi surface has perfect particle-hole 
symmetry, the vertex corrections are perfectly canceled 
out, and at approximately the half-filled state, owing to 
approximate particle-hole symmetry, vertex corrections 
are approximately canceled out and the values of Amax 
are suppressed. Here, the normal self-energy term corre- 
sponding to Fig. 2(c) make the mass enhancement fac- 
tor small. When the Fermi surface have perfect particle- 
hole symmetry, the terms corresponding to Figs. 2(b) and 
2(c) are perfectly canceled out. However, for n < 0.82, 
particle- hole symmetry deteriorates, and the mass en- 
hancement factor is much smaller than unity. Therefore, 
reliable numerical calculation cannot be obtained in the 
range of n < 0.82. 

In Fig. 7 we show the to dependence of Amax- With 
decreasing to, the spin-triplet state becomes dominant, 
and with increasing to, the spin-singlet state becomes 
dominant. If to is small, xo{q,0) have the character of 
one-dimensionality. In this case, the values of Amax are 
suppressed by the conflict of the peaks of xo(q,0) like 
the above case. 

In conclusion, we have investigated pairing symme- 
try and the transition temperature on the basis of a 
quasi-one-dimensional Hubbard model. We have solved 
the Eliashberg equation using the third-order perturba- 
tion theory with respect to the on-site repulsion U. We 
flnd that if n is shifted from the half-fllled state, the 
transitions into unconventional superconductivity is ex- 
pected. If one- dimensionality is weak, a spin-singlet pair- 
ing is more stable than a spin-triplet one. In contrast, if 
one-dimensionality is strong and n is far from the half- 
fllled state, a spin-triplet pairing is more stable than a 
spin-singlet one. Thus, we suggest the possibility of un- 
conventional superconductivity in /3-Nao.33V205 caused 
by the on-site Coulomb repulsion. 

Numerical calculation in this work was carried out at 
the Yukawa Institute Computer Facility. 
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